Transport through a disordered topological-metal strip 
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Features of a topological phase, and edge states in particular, may be obscured by overlapping 
in energy with a trivial conduction band. The topological nature of such a conductor, however, is 
revealed in its transport properties, especially in the presence of disorder. In this work, we explore 
the conductance behavior of such a system with disorder present, and contrast it with the quantized 
conductance in an ideal 2D topological insulator. Our analysis relies on numerics on a lattice system 
and analytics on a simple toy model. Interestingly, we find that as disorder is increased from zero, 
the edge conductivity initially falls from its quantized value; yet as disorder continues to increase, 
the conductivity recovers, and saturates at a value slightly below the quantized value of the clean 
system. We discuss how this effect can be understood from the tendency of the bulk states to 
localize, while the edge states remain delocalized. 



I. INTRODUCTION 

To date, topological behavior has been observed in 
many systems in zero magnetic field, usually by iden- 
tifying topologically induced edge states. Much of 
the theoretical attention has been given to topological 
insulators, 1 7 while many of the experimental observa- 
tions are of systems that are metallic, due to either dop- 
ing, or midgap statesP^It is natural then to ask - what 
kind of topological behavior can a metal, or a gapless 
system, exhibit? 

Several groups have studied the so-called topologi- 
cal conductor. Ref. [12] used a 2d Kane-Mele (or Hal- 
dane) topological insulator, with interstitial sites in each 
hexagon forming an additional trivial band to realize a 
topological conductor (see Fig. [I] for illustration). The 
'parasitic' metallic band was made to overlap with the 
topological-band edge states. Hybridization between the 
topological and parasitic metallic bands changed the 
edge-state spectrum in a peculiar way. Edge states 
which overlapped in both energy and parallel momen- 
tum with bulk metallic states not only did not disap- 
pear when hybridization was introduced, but rather, they 
doubled: Exact edge states emerged at energies above 
and below the metallic band, and a finite width spectral 
resonance remained where the unhybridized topological 
edge state used to be. An additional topological con- 
ductor system was studied by Bergman^ and Barkeshli 
and QiP^ They showed that a magneto-electric axion re- 
sponse (with time-reversal broken at the surface) per- 
sists even when the 3d topological insulator is doped, 
and Fermi-surfaces appear in its bulk, although it is no 
longer quantized. Similarly, Ref. 15 showed that a 3d 
topological metal still supports a special surface plasmon 
mode. 

The transport properties of a topological conductor, 
especially in the presence of disorder are the focus of this 
paper. These aspects of the topological conductor were so 



far mostly ignored. Several groups, however, have inves- 
tigated the effects of disorder on topological insulators. It 
was clearly shown that a sufficient amount of disorder will 
close a topological gap, resulting in a met alP^H Further- 
more, it was found that disorder could even induce topo- 
logical behavior in trivial semiconductors with spin-orbit 
coupling! 25 " 27 Only Ref.[28|so far also considered disorder 
effects in a topological-metal regime, and found that it 
does not qualitatively affect the non-universal magneto- 
electric effects found by Refs. I13|141 

Our study concentrates on the 'parasitic metal' fla- 
vor of the topological conductor which was described 
above Specifically, we investigate the conductance of 
a topological conductor strip as a function of disorder, 
with its chemical potential in the energy range where 
edge states and the parasitic metallic band coexist. In- 
creasing disorder will eventually localize the bulk states, 
and then we expect there to be rather little mixing be- 
tween the bulk states and the edge states, since the wave- 
function overlaps are exponentially small. In contrast, for 
weak disorder, where the bulk states are still mostly delo- 
calized, the wavefunction overlaps are much more signifi- 
cant, and we expect that no state will be distinguishable 
as an edge state. 

Our results, however, are quite surprising. The edge 
states appear to retain their significance for essentially 
all disorder strengths. Their presence and distinction 
is manifested in the persistence of a single highly con- 
ducting channel which survives until the disorder is suf- 
ficiently strong and the topological band associated with 
the honeycomb subsystem is destroyed. The distinct fea- 
ture of this effect is the appearance of a conductance 
minimum for the most conducting channel at some finite 
disorder (see Fig. |3|. We characterize this conductance 
minimum in systems of varying width and length, as well 
as for several values of the model parameters. 

The organization of the paper is as follows. In Sec. [TT] 
we describe in detail the parasitic metal model we 
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FIG. 1: A two dimensional lattice model of a topological 
conductor. The model is based on the Kane-Mele modeP^ 
on the honeycomb lattice, which is denoted by dashed lines. 
The Kane-Mele model consists of nearest-neighbor hopping t\ 
(green) between the two (A, B) sublattices of the honeycomb 
lattice (denoted by empty and filled circles respectively), and 
a complex second-nearest-neighbor hopping (a spin orbit cou- 
pling term) with opposite sign when clockwise (+iA, red) and 
counterclockwise (— iA, blue), as indicated by the curved ar- 
rows in the figure. In addition to sites of the honeycomb 
lattice, we include a new set of sites (C) at the centers of the 
hexagonal plaquettes of the honeycomb lattice, denoted by 
cross marks. The C-sites form a triangular sublattice, and 
nearest-neighbor hopping between them £2 (brown) forms a 
metallic band. To explore the interplay between the topo- 
logical insulators helical surface states and the bulk metallic 
band we mix the two systems by allowing hopping between 
the C-sites and the honeycomb lattice sites, £3 (purple). 



study. We describe the transfer-matrix method we use 
in Sec. |III| Our results, demonstrating the conductance 
minimum, as well as its dependence on the system param- 



eters are recounted in Sec. [IV| A qualitative understand- 
ing of the effect can be obtained using a simple model 
which we describe and analyze in Sec. [V| before making 
our concluding remarks in Sec. VI 



II. THE KANE-MELE PARASITIC BAND 
MODEL FOR TOPOLOGICAL CONDUCTOR 

We perform the disorder and transport analysis on a 
specific model for the topological conductor. We expect 
that the qualitative behavior will be independent of the 
specific model we employ. As mentioned in the introduc- 
tion, we will use the Haldane model with interstitial sites 
added. We will also consider the full Kane-Mele model 
with Rashba interaction included and find a similar be- 
havior as for the Haldane model. 

The Haldane modeP^ consists of spinless electrons 
on the honeycomb lattice, with nearest-neighbor hop- 
ping which gives the band structure of Graphene, and 
imaginary second-nearest-neighbor hopping which opens 
a gap, induces a non-zero integer Chern number in the 
conduction and valence bands, and produces chiral edge 
states. To this we add a parasitic metallic band, formed 
by an overlayed triangular lattice of sites, which occupy 



FIG. 2: (a) Energy spectrum of the topological conductor 
without disorder. Bands at the top and bottom are the bulk 
states of the honeycomb lattice. The band in the middle is 
the metallic bulk from the interstitial sites. The parameters 
in units of t\ are A = 0.1, £2 = 0.1, £3 = 0.1, and the width 
of the strip is M — 30 zig-zag lines, (b) Zoom of the band 
structure in (a). Different Fermi-energy regimes are indicated. 
The region between (II) and (III) has the same properties as 
region (I). 



the center of the honeycomb hexagons (see Fig. [T]). Hop- 
ping between the interstitial triangular lattice sites pro- 
duces a single topologically-trivial band, which overlaps 
in energy with the edge states of the Haldane model. We 
then allow the two subsystems to hybridize. The Hamil- 
tonian describing the combined model is 
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where z, j denote the sites of the composite honeycomb 
and interstitials lattice. The operators a, b denote the 
fermion annihilation operators on the two honeycomb 
sublattices (A and B), and c is the fermion annihila- 
tion operator on the triangular lattice sites C at the 
centers of the honeycomb plaquettes. The coefficients 
Uij = ±1 determine the sign of the imaginary second- 
nearest-neighbor hoppings, and are defined in Refs. 11121 
as Uij = -j=z • (djk x dki), where dkj is the unit vector 
pointing from site j to fc, and k is the intermediate site 
between j and i. The parameter t\ denotes the nearest- 
neighbor hopping on the honeycomb lattice, while £2 de- 
notes the nearest-neighbor hopping on the triangular lat- 
tice. Finally, £3 is the hybridization hopping between the 
honeycomb and triangular lattice sites. The random on- 
site potential Vi on every site of the lattice (honeycomb 
and triangular lattice sites) has a uniform distribution 
between [-W, W}. 

The single spin system we analyzed consisted of zig-zag 
strips of the topological conductor with the Hamiltonian 
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Q. Its generic band structure is shown in Fig. |2j In- 
deed we see that the metallic band intervenes between 
the conduction and valence bands of the Haldane por- 
tion of the band structure. Indeed we see regions in the 
spectrum where in a single value of the momentum paral- 
lel to the edge, there are two edge states at energy above 
and below the metallic band. This is a manifestation of 
the exiling effect pointed out in Ref. 12, and indicates 
that the metallic band appeared on top of the edge-states 
branch in this region, and hybridization expelled the edge 
states while doubling them. 



A. Rashba coupling 

We have also studied the effect of adding a Rashba 
type spin-orbit coupling on the honeycomb lattice to 
our model 0, with the electron spin— \ restored. The 
Rashba interaction mixes between the two spin flavors as 
follows 

U R = l\ R ^ a \a^oc(3 X dij) z bjp + h.c, (2) 

(ij)a/3 

where s is the vector of Pauli matrices for the electron 
spin, dij is the unit vector pointing from site j to z, and 
a, /3 are spin indices. 

Even with Rashba interaction included, we find the 
same qualitative behavior as in the spinless case de- 
scribed above. 



B. Fermi-energy regimes 

The systems' transport properties depend closely on its 
Fermi energy. The range of Fermi energies between the 
conduction and valence bands of the topological subsys- 
tem can be split into three important regions. Region I 
denotes the case where the edge states and metallic bands 
overlap in energy but not in momentum, and therefore 
coexist in the clean limit. Region II denotes the energy 
range of the parasitic metal where no edge states appear. 
Region III is the energy range that has only edge states. 
These regions are indicated in Fig. [2^b) for the spinless 
model. 



current flowing through the system in a two-terminal 
device p3 

G=^|Ti»| 2 , (3) 

n,m 

where T^(Ep) is the transmission coefficient for going 
from mode m in the left lead to mode n in the right 
lead at the Fermi energy, Ep. To ensure probability cur- 
rent conservation, it is important to note that the current 
associated with a scattered wave is proportional to the 
square of the wavefunction multiplied by the velocity. 
Therefore T nm is given by the outgoing current ampli- 
tude, i.e., the wave amplitude times the square root of 
the velocity of an electron leaving the device in mode n 
through one lead if the incoming current amplitude in 
mode m in the other lead is set to unity. 

Eq. (§ can also be written as the sum of conductances 
per channel g m : 

G=^^> m , (4) 

m 

where g m = \T^\ 2 is the conductance of channel 
m, given by the probability that an electron entering 
the system in mode m is transmitted through the con- 
ductor. We find the transmission coefficients T^(Ep) 
using the transfer-matrix method, as described in Ap- 
pendix [] We note that this method, particularly when 
disorder is concerned, has an instability which restricts 
the system size possible. This instability arises due to 
the existence of imaginary momentum modes at the de- 
sired energy range which appear because we are exploring 
the conductance of the topological conductor for energies 
in the bulk gap of the honeycomb lattice. Indeed, the 
number of imaginary momentum modes due to the hon- 
eycomb subsystem increases linearly with the number of 
zig-zag lines in the strip, thus restricting the accessible 
system size!^ Nevertheless, the transfer-matrix method 
provides reliable results for a range of strip sizes. These 
results could, in principle, be improved upon by using 
S-matrix methods. 34 



III. LANDAUER FORMALISM FOR THE STRIP 

In order to analyze the two-terminal transport through 
the topological conductor, we use the Landauer-Biittiker 
formalism! 30 " 32 1 We envision our system as consisting of 
a long strip of the topological conductor, which is dis- 
ordered in a finite region (Fig. 



14). The regions to the 



left and right of the disordered region form the ballis- 
tic leads, with several transverse modes. The Landauer- 
Biittiker formalism gives the following formula for the 



IV. RESULTS 

We investigated the conductance of the topological 
conductor strip for different Fermi-energy regimes and 
analyzed its dependence on system and model parame- 
ters. For Fermi energies where metallic bulk and edge 
states coexist (region I), we find a distinct minimum in 
the conductance of the most conducting channel as a 
function of disorder strength, followed by a revival to- 
wards ballistic transport in this channel. This is a sur- 
prising feature given the mixing between bulk and edge. 
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FIG. 3: Illustration of the 'dip and recovery' feature of the 
edge channel conductance. Conductance g n of different chan- 
nels is plotted as a function of disorder strength W . The 
edge channel conductance shows a sharp minimum g m i n at 
finite disorder Wo and a broad maximum g ma x- Parameters 
are E F 0.1, A = t 3 = 0.1, and M — L — 20. 



A. Region I 

The conductance for each channel for a Fermi energy 
in region I is shown in Fig. [3j The generic feature we find 
in this regime is a single highly conducting channel that 
persists up to large disorder and decays when disorder is 
presumably strong enough to close the gap of the honey- 
comb subsystem. For small disorder there is a sharp dip 
in the conductance followed by a broad maximum. 



1. Dependence on system size 

In general, increasing the width of our strip has two 
opposing effects. One, the overlap of edge states on op- 
posite sides of the strip is reduced and thereby the con- 
ductance should increase. Two, a wider strip hosts more 
bulk channels and leakage from the edge into the bulk 
is increased which reduces the conductance. Increasing 
the length of our strip increases the probability of elec- 
trons inside the wire being backscattered or leaking into 
the bulk since they travel a longer distance inside the 
disordered region. This should decrease the conductance 
while the net effect of changing the system width is not 
apparent. 

The dependencies of the characteristic values of the 
conductance, the minimum g m i n , the maximum g max at 
higher disorder, and the disorder strength at the min- 
imum Wo (cf. Fig [3]), on the system size are given in 
Fig. [4j We find that these quantities have rather differ- 
ent dependence on the system dimensions. g m [ n shows 
an exponential decay with increasing system length and 
a small increase with increasing width. The disorder 
strength Wo at the minimum seems to be independent 
of the system width but shifts to smaller disorder val- 
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FIG. 4: Dependence of the conductance dip and maximum on 
system size (Region I), (a)-(c) Dependence on length L of the 
system for M — 20. Note the log-scale on the y-axis in (a); the 
dashed line is a guide to the eye to emphasize the exponential 
dependence, (d)-(f) Dependence on system width for (d), (e) 
L = 40 and (f) L 20. Parameters: E F 0.1, A = t 3 = 0.1. 
While the observed feature is nearly independent of the sys- 
tem width, the conductance minimum decreases exponentially 
with system length. The maximum conductance is nearly in- 
dependent of the system length and the deviation from the 
quantized value presumably results from backscattering at the 
interface between leads and disordered region. 



ues with increasing system length. g m i n , #max, and Wo 
are (nearly) independent of the strip width which shows 
that either the two opposing effects of changing the width 
are not very significant or that they (nearly) cancel each 
other for the studied parameter region. 

Rather importantly, it appears that the most conduct- 
ing channel recovers at high disorder and becomes effec- 
tively ballistic. The conductance maximum g max seems 
to be independent of the system length. Its deviation 
from the quantized value is presumably the result of 
backscattering at the interface between leads and disor- 
dered wire rather than from effects inside the disordered 
region. 



2. Dependence on system parameters 

An increase in the coupling £3 between the honeycomb 
subsystem and the triangular lattice increases the prob- 
ability of electrons leaking from the edge into the local- 
ized metallic bulk states which should reduce the con- 
ductance of the chiral edge states. Imaginary second- 
nearest-neighbor hopping A is the essential ingredient for 
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the hopping parameters A and £3 is illustrated in Fig.[5|a) 
and (b). We find that the characteristic shape of the 
conductance is left unaffected by changes in the hopping 
parameters of over a wide range of values. As expected, 
an increase in the hybridization £3 decreases the conduc- 
tance while an increase in the second-nearest-neighbor 
hopping A increases it and leads to a broader maximum 
due to the larger bulk gap of the honeycomb subsystem. 
One might think that the dependence of the conductance 
maximum on the hybridization £3 contradicts our earlier 
conclusion that the deviation of the conductance maxi- 
mum from the quantized value results from backscatter- 
ing effects at the interface between conductor and leads. 
This is, however, not the case since an increasing coupling 
between the triangular and honeycomb lattices increases 
backscattering everywhere, including at the interface be- 
tween conductor and leads, leading to the observed be- 
havior. 

When varying the Fermi energy within region I to val- 
ues away from the metallic band (see Fig. [5jc)), first the 
conductance of the most conducting channel increases 
over the entire disorder range because of a decreasing 
number of metallic states but the characteristic shape is 
unaffected. However, for Fermi energies too close to the 
bulk of the honeycomb lattice, the conductance starts to 
lose its shape. Even small disorder can in this case de- 
stroy to topological bands through partial closing of the 
bulk gap of the honeycomb subsystem. 




B. Region II 
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FIG. 5: Region I. Dependence of the edge channel conduc- 
tance on system parameters. Dependence on (a) hybridiza- 
tion £3 (A = 0.1, M = 20) and (b) second-nearest-neighbor 
hopping A (£3 = 0.1, M = 10) for Ef — 0.1 and L = 20. (c) 
Dependence on Fermi energy with A = £3 = 0.1, M — 10, 
and L — 20; (d) indicates the energy range from E = 0.1 to 
E — —0.3. Decreasing the hybridization £3 and increasing the 
second-nearest-neighbor hopping A increases the conductance 
over the entire disorder range. Moving the Fermi energy away 
from the metallic bulk has the same effect until i^F is too close 
the honeycomb valence band when the characteristic feature 
is starting to vanish. 



a topological phase as it opens a gap and produces chi- 
ral edge states. An increase in A widens the gap and 
thereby makes the system more robust against disorder 1 
which should increase the conductance and broaden the 
maximum. 

The dependence of the edge channel conductance on 



In this region, the Fermi energy lies deep within the 
(trivial) metallic band. Without hybridization between 
the honeycomb and the triangular lattice, the edge states 
would traverse the bulk gap unhindered through the 
metallic band. For finite hybridization, however, the edge 
states are expelled out of the metallic band region to 
higher and lower energies (cf. Fig. [2J. We would there- 
fore expect the conductance in this regime to show bulk 
state behavior, i.e., rapid decay with increasing disorder. 
Decreasing the hybridization should recover the charac- 
teristic edge state behavior. We find, however, that even 
for strong hybridization £3 = 0.1, for a certain range of 
parameters there is a conducting channel which shows 
the characteristic behavior of the edge channel (Fig. [6]) . 
This highly conducting channel exists only for a range of 
Fermi energies within the metallic band. 

The explanation of the highly conducting channel in 
this regime is consistent with the notion of leftover spec- 
tral resonances of the edge states, once they are absorbed 
by the bulk. When the surface bands are pushed away 
from the metallic band they leave behind a finite lifetime 
surface resonance in place of the original surface states 
(these are referred to as 'ghosts' in Ref . H2j) . 



FIG. 6: Region II. Conductance of most conducting channel 
for a Fermi energy deep within the metallic band, (a) Depen- 
dence on disorder strength for different hopping strengths £3. 
E F = 0.2, A = 0.1, M = 20, and L = 20. (b) Conductance of 
each channel at different Fermi energies. For Fermi energies 
between E = 0.21 and E = 0.22 a single highly conducting 
channel exists possibly connected to surface resonances. Pa- 
rameters are W 1.3, A = £3 = 0.1, M = 10 and L = 20. 
(c) Conductance as a function of disorder at Ef = 0.21 for 
A = £3 = 0.1, M = 10 and L = 20. A single channel stands 
out, highly conducting at finite disorder. Dashed line indi- 
cates the disorder strength taken in (b) . (c) shows the energy 
range from E = 0.18 to E = 0.23. 



C. Region III 

In this region the Fermi energy lies above the metallic 
band and below the conduction band of the honeycomb 
subsystem. The conductance is only determined by the 
edge states and we would expect a quantized value up to 
a certain critical disorder strength as in the case without 
a metallic band. However, for Fermi energies very close 
to the metallic band there is a residual effect from the 
metallic bulk states and the conductance shows a small 
dip in the otherwise quantized conductance (Fig. [7]). 

Disorder closes the bulk gap of the honeycomb subsys- 
tem by creating new bulk states above the valence and be- 
low the conduction band, and the conductance decreases 
rapidly with disorder for larger Fermi energies. In addi- 
tion, disorder also creates new metallic states above the 
metallic band. These new states then coexist with edge 
states as in region I which explains the observed small 
dip at small disorder. 



D. Rashba coupling (Region I) 



FIG. 7: Region III. Conductance of the edge state for different 
Fermi energies. Parameters are A = £3 = 0.1 and M — L — 
20. Inlet indicates the energy range form E — 0.35 to E — 
0.45. At Fermi energies close to the metallic bulk (Ef — 0.35), 
disorder induced states above the bulk band still lead to a dip 
in the conductance. The closing of the honeycomb bulk gap 
due to disorder quickly decreases the conductance in the entire 
energy range. 




FIG. 8: Band structure with Rashba coupling (a) Xr = 0.1 
and (b) Xr — 0.3. Parameters are A — £3 — 0.1, and M — 
30. Rashba coupling reduces the band gap of the honeycomb 
subsystem. 



small Rashba coupling does not change the band struc- 
ture dramatically. For large Rashba coupling the edge 
bands are clearly separated in momentum and the struc- 
ture of the metallic band changes significantly. The bulk 
gap of the honeycomb subsystem is reduced such that re- 
gion III vanishes. However, we find that Rashba coupling 
does not affect the characteristic shape of the edge chan- 
nel conductance as seen in Fig.[9| Increasing Rashba cou- 
pling decreases the conductance over the entire disorder 
range but leaves the distinct 'dip and recovery' behav- 
ior of the conductance intact. As for the model without 
Rashba coupling, decreasing the hybridization £3 between 
honeycomb and triangular lattice increases the conduc- 
tance. 



We also considered a spinful version of ([I]) with Rashba 
type spin-orbit coupling included between nearest- 
neighbor sites of the honeycomb subsystem We con- 
sidered the effects of spin-independent disorder only. As 
shown in Fig. |§J aside from the doubling of the bands, 



V. INTERPRETATION AND TOY MODEL 

The nature of the conductance dip found in the previ- 
ous section can be qualitatively understood using a sim- 



7 




FIG. 9: Region I with Rashba coupling. Conductance of one 
edge channel for (a) £3 = 0.05 and (b) £3 = 0.1 for different 
Rashba coupling strengths Xr. Parameters are Ep — 0.05, 
A = 0.1, M — 10, and L — 20. We plotted only a single spin 
orientation for readability. Rashba coupling overall reduces 
the conductance. The effect is enhanced for larger hybridiza- 
tion between the honeycomb and metallic subsystems. 



pie toy model. Our toy model consists of a single chiral 
mode that is mixed, due to disorder, with multiple bulk 
modes. This model neglects inter-edge scattering. We 
further approximate the most conducting channel from 
one side to the other as the sum of two incoherent contri- 
butions: edge-to-edge g e -e — \ \T e ,e\ 2 (through just the 

edge channel) and bulk-to-edge g e - b = ^ Y^m^e \ T e,™\ 2 
transport. 

The transmission coefficient between mode m in the 
left lead to mode n in the right lead, can be computed 
from the retarded Green's function connecting these two 
modes 33 



ihy/v 



(5) 



where v n is the electron velocity in mode n at the Fermi 
energy. For the sake of brevity we will take H = 1 from 
this point forward. 

Indeed, since this system has a chiral mode, strictly 
speaking, it will always have a total conductance bigger 
than the single channel, e 2 /h. Nevertheless, we will only 
be interested in the conductance involving the bare chiral 
edge states. So if the chiral mode manages to wander 
deep into the bulk due to the disorder, we will ignore its 
conductance. This is consistent with the numerics, since 
under these circumstances one would no longer be able 
to separate the system into two independent parts, each 
involving a single chiral mode. 



In the following we study the toy model, explicitly cal- 
culate the conductance for different system lengths and 
widths, and use the results to explain the origin of the 
conductance dip. 



A. Toy model and Green's functions 

Consider a system with a single chiral edge mode and 
multiple bulk modes. We assume that the chiral edge 
mode can only hybridize with the bulk when impurities 
are introduced. The hybridization strength will serve as 
a proxy for disorder strength. 

We describe our system with the action 



duo 



dxdx' 



X \x)G- 1 (x-x')x(x') 



I 2k 

+ i>i(x)F-^(x - x')ip m (x') 
- \xHx)V n (x)ip n {x) + h.c] 



(6) 



where the field x( x ) represents an electron at position 
x in the bare edge mode, and the fields ijj n (x) represent 
an electron at x in bulk mode n. In addition, G and F 
denote the bare edge and bulk state propagators, respec- 
tively, and V denotes the random potential scattering 
between the edge state and bulk channels. We will adopt 
the uniform distribution of random potential strength. 
The bare propagators are translationally invariant in the 
direction through the length of the strip (hence the de- 
pendence only on the relative coordinate x—x'\ while the 
impurity scattering V(x) is completely local. All terms 
depend on frequency cj, which for the sake of brevity has 
been kept implicit. The Einstein summation convention 
has been used for the bulk channel indices n, m. In mo- 
mentum space the form of the propagators is 



G 1 (k,oj) = uj + if] — v e k 



u> + ir] 



2M r 



(7) 



where 77 > is infinitesimal and positive, and A4 n and 
/i n are the effective mass and the chemical potential, re- 
spectively, in bulk mode n. Fourier transforming from 
momentum space to real space we find 



/-Too 
G(k,E F )e ik: 
-CO 

/ + OO 
F n , m (k,E F ) 
-00 



x 0(x) 



(8) 



_ p ikb,n\x\£ 

u n,m 1 



where k^^ = ^2j\A n {E F + fi n ) is the Fermi momentum 
in the bulk mode n, and k e = The chiral mode prop- 
agator G reflects a single propagation direction. G(x) is 
zero when x < 0, and this reverses when we take 77 < 
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for the advanced propagator. For the bulk modes, we 
assume F is diagonal, neglecting the mixing of different 
bulk conduction channels. 

The quantity V(x) is not the bare random potential, 
but rather matrix elements between bulk modes and the 
edge state, due to disorder scattering, and we therefore 
assume they are uncorrelated random complex variables. 
The random phase captures the random location of the 
focal points of impurity induced edge-bulk scattering, 
and makes sure our results are not affected by momen- 
tum conservation. In this manner we take into account 
the diffusive nature of the bulk, while still using the trans- 
lationally invariant form of the bulk propagator F. 

We discretize the x coordinate into i = 1 ... N, as il- 
lustrated in Fig. [lO} and set V n (i = 1) = V n (i = N) = 
to represent clean leads so that we can still separate in 
our calculations the edge conduction channel from all the 
others. For i ^ 1,7V we take V n (j) = tj^ n e^^ n . Now 
G and F are both square matrices with sizes N x N 
and NM x 7VM, respectively, where M is the number 
of channels (and the width of the lattice model of our 
strip). In contrast, V n is a non-square matrix of dimen- 
sions N x NM. Discretizing x requires one subtle change 
to the chiral propagator G in ^ 



G(xi — Xj, Ep) = - 



I 



ik b) \xi-Xj\? 

u n,m 5 



(9) 



where the step function 6{x) at x = is defined to be 
zero. To find the transmission coefficients we seek, we 
will need to calculate the Green's function for electrons 
starting out at the left lead i = 1 and coming out at the 
right lead i = TV, at the Fermi energy 

Ge,e( X = N,X* = 1) = (Xi=NXi=l) 
Ge,m(x = N,x' = 1) = (Xi=N^m,i=l) • 

where we specifed two contributions: electron from edge 
to edge mode, Q e ^ ei and electron from bulk mode m to 
edge mode, G e ,m- With the position index suppressed for 
the sake of brevity, our action is of the form 



(10) 



duj 



G' 1 -V 
_yt f-1 



(11) 



Next we want to invert the matrix, to find the compo- 



nents of the full Green's function Q 
We find 

1 



G' 1 -V 
-W F -i 



e ' e G- 1 - VFV^ 

0m,e Gm,nV n G 
1 

F- 1 - V^GV 



(12) 



where we have used the Einstein summation convention. 
The non-square matrices G e ,m and Q m ,e are the Green's 
function elements that capture the bulk-edge and edge- 
bulk mixing, respectively; Q e ^ e and Q n ,m are the renor- 
malized Green's functions for the edge and bulk modes, 
respectively. 

As we will see in what follows, the simple toy model 
already affords us a qualitative understanding of the con- 
ductance dip, and its dependence on the system dimen- 
sions. 

B. Conductance dip of different system size 



i=l 


=2 i=3 




i=N-l 


i=N 






jr. 
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... V V 


^ Bulk model 























FIG. 10: Illustration of the toy model. Bulk modes are inde- 
pendent of each other since we assume that the bulk channels 
describe the eigenmodes of the diffusive bulk. The edge mode 
couples to the bulk through (N — 2) scattering sites, giving in 
total (N — 2)M couplings V n (i). The coupling strength and 
phase are uniformly distributed. 



1. Single contact with single bulk case 

It is instructive to first consider the case of a single 
bulk mode m = b and a single impurity. In this case 
we take N = 3 and M = 1, and for simp licity we take 
k e = kb = kp-> and V — te 1 ^. From Eqs. ( 9|l2 ), we can 
show that 



e ik F L 


1 


v b 1 + 


t' 2 /2 


e ik F L 1 _ 


t' 2 /2 


V e 1 + 


-t' 2 /2 


e ik F L 


t' 


^V e V b 1 - 


f- 1' 2 /2 


e ik F L 


a 



(13) 



where L = (N — l)a = 2a, a is the lattice spacing, and 
t' = t/^v e Vb is a dimensionless coupling strength. From 
these we find the conductances 

in 2 



e 
~h 

h 



2 r 



l-t' 2 /2 



1 + 1' 2 /2 

a 

l + t' 2 /2 



(14) 
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We can now easily see that g e -e — \ for both the clean 
limit t' = and the infinite disorder limit t' — >• oo. In 
between there is a minimum of g e -e — at t' = y/2. In 
contrast, g e _5 starts at zero for t' = 0, and rises to a 
maximum value of g e -6 = at t' = a/2, before decay- 
ing back down to zero. Note that the Fermi velocities 
only rescale the disorder strength, and did not otherwise 
influence the conductance behavior. 



where V n (xj) = e ik * x iV n (xj) /2, such that 

(VFV^ij = En, m Vn(i)Fn,m(hj)Vm(jy. The deriva- 
tion of this result is somewhat involved, and so we defer 
its details to a future publication.^ More generally, for 
j ^ 1, we can show that 

g (Nj )-itli e ^ N -j) ^-VF^) jth 




where (A)jth indicates the j row of matrix is re- 



placed by a vector (a) 



FIG. 11: Conductance contribution via different channels as 
a function of coupling strength (x-axis) for the system of a 
single chiral and bulk mode with a single coupling: from edge 
to edge g e -e, edge to bulk gb-e, bulk to edge g e -b, and bulk to 
bulk gb-b- The edge to edge conductance has a minimum at 
finite disorder before recovering to the quantized value. The 

2 

conductance is in units of 4r- 



Finally using Eq. (|| we find all the different con- 
tributions to the conductance (edge/bulk-edge/bulk), as 
Fig. pTj^a) shows the total conduc- 
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plotted in Fig. 

tance of the strip, consisting of adding up the edge-to- 
edge and bulk-to-edge channels. Fig. [Tljb) shows the 
possible individual conductances. Two points are note- 
worthy: (i) In the strong disorder limit, the conductance 
of edge-to-edge channel converges to unity, (ii) At the 
conductance dip minimum, the biggest contribution to 
the conductance comes from the bulk-edge channel. In 
other words, at strong disorder, the edge-edge channel 
conductance recovers to unity even without accounting 
for the localization of the bulk, and the size dependence 
of the conductance near the dip will arise from the bulk- 
edge channel contribution, and will be rather indepen- 
dent of the edge-edge channel behavior. 



2. Multiple bulk modes and impurities 

Consider now a general system of length TV (i.e., hav- 
ing TV — 2 scattering sites), and width (bulk mode num- 
ber) M. The Green's function components we need to 
calculate to find the transmission coefficient G{N, 1) can 
be found in closed form. In particular, the edge-edge 
Green's function is found to be 



^(iV,l) = -^^l) 



det(G T - VFVt) 
det(G + VFW) 



is 



(0, — i, — i, . . . , —i, 0). Concomitantly, the Green's func- 
tion bulk-edge mixing component is 

N-l 

<5 e , m (7V, 1) = J2 E M N > j)Vnti)Fn,mtii 1) (17) 
n j=2 

which is the sum of the (N — 2) Green's functions from 
clean bulk to renormalized chiral mode. 



3. Length dependence 

The conductance minimum g m i n is expected to drop 
exponentially with the system's length, as shown in Fig. 
|4ja). The toy model's results are shown in Fig 12 
Within our toy model, the length dependence of g ] 
attributed to the mixed edge-bulk contribution to the 
conductance g e -b- Indeed, we approximate the edge- 
related conduction as an incoherent sum of two contribu- 
tions: conductance of electrons entering and leaving the 
strip in an edge state, g e -e: and entering as bulk states, 
but leaving the strip as edge electrons, g e -b- Even with 
multiple scatterers, g e _ e behaves as Eq. (15) indicates: it 
drops down from unity (in units of ^-) to zero at moder- 
ate disorder, then rises back to unity at strong disorder. 
In Fig. [l2]fo) we see that the minimum g e - e is indepen- 
dent of the system length. Therefore, the only possible 
source of the g m { n dependence on the system length must 
come from g e -b, which starts at zero for the clean limit, 
rises for moderate disorder, then peaks and decays as we 
go to strong disorder (see Fig. [l2^c)). The inset in Fig. 
[l2^ a) shows the conductance # e -& depends exponentially 
on the number of scattering points, and hence on the 
length of the system. 



(15) 



4- Width dependence 

We now explore the width dependence of the conduc- 
tance through the edge states. To this end, we can use 
a simple approximation VFV^ ~ (VFV^). This stems 
from the fact that ( (VFV^) n ) / ((VFV^)) n -> 1 when 
the number of bulk modes goes to infinity M — >> 00. This 
holds as long as the distribution of the random matrix 
elements Vj yU has a finite variance. Our approximation is 
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5 10 15 20 25 5 10 15 20 25 

w w 

FIG. 12: (a) Sum of the conductances of chiral-to-chiral (b) 
and bulk-to-chiral (c) channels for M=4 and N=5 (blue), 7 
(green), 9 (red) and 11 (cyan) for k e = anc ^ ^ b ~ ^T"- 
The inset in (a) shows the exponential decay of the minimum 
conductance with system length in agreement with the nu- 
merical calculation in Fig. H[a). The conductance is in units 





w w Vm 



FIG. 13: (a) Sum of conductance of chiral-to-chiral (b) 
and bulk-to-chiral (c) channels for N=5 and M=4 (blue), 7 
(green), 10 (red) for k e = ±£l an d k b = ^l. The conduc- 
tance minimum is unaffected by changes is the system width 
in agreement with the results of the numerical calculation, c) 
Renormalizing the disorder strength by \/M, g e -b, which is 
the dominant contribution to the conductance in the region 
of the minimum, becomes independent of the system width. 

2 

The conductance is in units of 



therefore very good when there is a large number of bulk 
modes M ^> 1, and we can then write 

(G (N i))~-^-i) ^ r -( Vfvt )) (18) 
{ y eA lM,l))- ie det(G + (VFVi)) ' ( ' 

A similar expression holds for Q e ^ e {N,j ^ 1) as well. 

The width dependence of the conductance can now be 
deduced from 

«VFVt» y . 

= (_1)*+VM^) J2K,m(iJ) (Vn(i)V m (jY) 
nm 

= ( _ l) i + i^e(x i -, j )^ F ^ (iii)( J.. C ^ (19) 
nm 

= %^ 2 5>„,„(M) ~ MSijW 2 , 

n 

where the variance of the disorder is the disorder strength 
W squared. From this we see that the dependence on M 
can be absorbed into W, 

W eff = ^fMW . (20) 

As a consequence we expect that increasing the width 
(increasing M) only shifts the minimum to different dis- 
order strength Wo, and does not change g m in- 

This simple dependence on M turns out to work quite 
well, when applied to the toy model. Fig. [l3[a) shows the 



conductance through the chiral mode for different system 
widths, M. The minimum conductances in the plots are 
indeed all similar, consistent with the numerical results 
shown in Fig. Qd). The conductance is dominated by 
the bulk-to-chiral channel (Fig. 13 'b)). Furthermore, 



scaling the disorder strength by vM, the three lines of 
Fig. fl3Fb) collapse onto a single line in Fig. 13 'c). Only 



the width of the minimum, and not its depth, is affected 
by the channel number. 



C. Summary of toy model analysis 

We have shown in this section that coupling between 
edge and bulk modes leads in general to the conductance 
dip behavior with increasing disorder strength. The de- 
pendence of the minimum conductance at the dip g m i n 
on the system size is consistent with the previous calcula- 
tion: it exponentially decays with the system length and 
is roughly independent of the system width. Neverthe- 
less, the toy model is limited. First of all, the localization 
of the whole system at strong disorder is not captured by 
it because the toy model assumes a fixed number of chiral 
and bulk modes at a certain Fermi level. Second, the toy 
model is directly concerned with the tunneling between 
modes and not equipped to capture reflection within a 
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single mode. As a consequence, the non-quantized value 
of the conductance plateau at strong disorder cannot be 
captured. Third, the toy model is not suitable to describe 
the region II calculation where edge and bulk modes are 
coupled already in the clean limit, and as such are not 
well separated contrary to what we assume in the toy 
model. 



VI. CONCLUSIONS 

In this manuscript we explored the transport signa- 
tures in a topological conductor phase. The particular 
realization of a 2d topological conductor we considered 
here, had an edge state originating in the high-energy 
bands, overlapping in energy (but not in momentum) 
with a trivial band at the Fermi energy (see the spec- 
trum in Fig. [2|. This model should suffice to capture the 
generic qualitative features of the non-protected topolog- 
ical aspects of topological conductors. 

To find topological signatures, we calculated the two- 
terminal conductance through a topological conductor 
strip in the presence of disorder in the form of a ran- 
dom on-site chemical potential. Our results clearly show 
that even when a bulk gap is absent, the edge states, 
which are no longer protected, maintain their individu- 
ality and dominate the transport in the strip. Inspecting 
the effects of disorder on all the conducting channels of 
the topological conductor, a single channel clearly stands 
out (see Fig. [3|. The most conducting channel in the 
strip has a conductance which initially decays rapidly as 
we increase disorder strength, but then reaches a mini- 
mum, and then recovers to a value consistent with ballis- 
tic transport through the disordered region. The dip and 
recovery features single out the topological edge chan- 
nel, and show that even when the system is gapless, its 
capacity for nearly dissipationless transport is restored 
once disorder suppresses transport through the bulk. 

The origin of the 'dip and recovery' features (Fig. [3| 
presumably has to do with localization effects in the 
bulk. Intuitively, as the bulk states become less trans- 
parent, the edge state, which cannot localize since it is 
chiral, reemerges to dominate transport. To gain a bet- 
ter qualitative understanding of the main features in our 
numerics, we developed a simple toy model, consisting of 
a single chiral mode that can randomly scatter to sev- 
eral regular channels. The model assumes that the bulk 
is diffusive, and does not invoke localization explicitly. 
Disorder strength is encoded in the random scattering 
strength between the chiral edge and the bulk. The dif- 
fusive nature of the bulk is taken into account by us- 
ing random phases in the edge-bulk scattering elements. 
Despite its simplifications, this model demonstrates that 
generically the edge-edge scattering contribution to the 
conducta nce g e - e shows the 'dip and recovery' behavior 
(see Figs. [TT|l2l ). 

Thus it appears the full-fledged Anderson localization, 
which is necessary to eliminate conductivity in the bulk, 



is not necessary for the revival of the edge state conduc- 
tance. Rather, the dip and recovery feature appears to 
be guaranteed once a separate chiral mode is present, as- 
suming it is unable to scatter to its counter-propagating 
partner. In the numerics, we do have both edge chan- 
nels, but because their wave functions essentially do not 
overlap, disorder can not produce direct scattering be- 
tween them. Note, however, that it appears that the 'dip 
and recovery' feature may not occur when the edge state 
overlaps in energy, and therefore, can hybridize, with the 
bulk states of its progenitor topological band. 

The unique transport features in topological conduc- 
tor strips will have additional signatures. First, we 
note that the results we presented here, when applied 
to a spinful electron system, rely on time-reversal invari- 
ance, and require that no scattering occurs between the 
counter-propagating, opposite-spin modes on the same 
edge. Therefore, when moderate disorder is present, and 
the strip conductance is dominated by the descendant of 
the edge channel, the conductance of the strip would be 
strongly suppressed by a magnetic (Zeeman) field nor- 
mal to the Rashba spin direction. Such a Zeeman field 
opens up a gap in the edge states, pushing them away 
from the Fermi energy, in which case the edge state no 
longer contributes to transport. The highly conducting 
mode we found would then no longer contribute to the 
overall conductance through the strip, which will now be 
purely diffusive, with no remnant of ballistic transport. 
Magnetic impurities at the edge would have a similar de- 
structive effect. Second, probing the current density in 
the sample would disclose that in moderate disorder the 
conduction is dominated by regions near the edge. Prob- 
ing these effects experimentally would rely on materials 
or quantum wells having a band structure similar to the 
one considered here: A metallic band coexisting at the 
same energy with the edge states of a different band. 

The 'dip and recovery' feature arising in transport due 
to disorder was shown here to appear in a non-interacting 
2d system. It is interesting to ask how these features 
would change when considering Coulomb interactions. 
This could also be explored within the toy model we pre- 
sented, by promoting the electronic modes to Luttinger 
liquids with interaction between the various modes. We 
could then also consider a more complete model contain- 
ing the chiral modes on both sides of the strip, and follow 
the mixing between them explicitly. Similarly, one could 
ask about similar resurgence of edge physics in a 3d topo- 
logical conductor. Instead of conductance, however, the 
appropriate topological signatures in 3d should be the 
Kerr effect rotation angle, which encodes the unique re- 
sponse of topological edge states to an oscillating electro- 
magnetic field (see also Ref . 123]) . 

The work presented in this manuscript adds to the no- 
tion that topological behavior has distinct and strong 
signatures even in the absence of protection due to 
a band gap. This observation, as articulated also in 
Refs. I13|14|21| significantly extends the class of systems 
(and materials) in which topological effects could be 
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found. 
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Appendix: Transfer-matrix method 

In this section, we will demonstrate a way to calculate 
the transmission coefficients of our system using the 
transfer-matrix method!^ 

Naively, all we have to do to calculate the transmission 
coefficients through our system is first find the basis of 
modes in the clean limit (the modes in the leads) ip n (x), 
then starting with an electron in just one of the (clean- 
limit) modes in the left lead ip n (x = au) (cf. Fig. [14] ), 
propagate the wavefunction through the dirty sample, 
and calculate its weight in the various modes of the 
(clean) right lead tj){x aw) Ysm T nm^rn(x aw). 
These weights are precisely the transmission coefficients. 

In practice, we run into some difficulties. Namely, the 
calculated ip(x = aw) contains weight in some unphys- 
ical exponentially growing modes ~ e +N^ not just in 
the propagating modes. As a consequence, we have to 
perform a more complicated numerical procedure. 

The main steps of the calculation are the following: 
After determining the transfer matrix, we calculate the 
eigenvalues and -vectors of the transfer matrix which cor- 
respond to the possible modes and wavefunctions in the 
leads. For computational convenience we start the calcu- 
lation in lead R. Propagating one of these wavefunctions 
backwards through the system using the transfer matrix 
gives us a wavefunction in lead L which is a superposition 
of all possible wavefunctions in the leads, i.e., a superpo- 
sition of all eigenvectors of the transfer matrix including 
unphysical (exponentially increasing) modes. Next, in 
order to determine the transmission coefficients, we need 
to find a superposition of only physical modes in lead 
L which, after propagation through the system, leads to 
a superposition of one right-moving mode with weight 
one and decaying modes in lead R (cf. Fig. 15). To 



find the necessary coefficients of the modes in lead L, we 
propagate a superposition of one right-moving mode with 
weight one and all decaying modes backwards through 
the system and choose the coefficients of the decaying 
modes such that the unphysical modes in lead L are 
canceled by destructive interference. Doing this for all 
right-moving modes gives us information about all trans- 
mission channels and we can determine the transmission 
coefficients for arbitrary modes n and m using prob- 
ability current conservation. 

The lattice model of the topological conductor is illus- 



trated in Fig. 14 The system is a honeycomb strip with 
a width of M zig-zag lines with an overlayed triangular 
lattice with sites that occupy the centers of the honey- 
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FIG. 14: Lattice system with a width of M zig-zag lines. 
The disordered region of length L is the conductor which is 
connected to a lead on either side. Dashed rectangles show 
the unit cells of the strip geometry. 



comb hexagons. The unit cell is chosen as explained in 



Fig. [14] with dashed rectangles indicating the unit cells. 
A region with a length of L unit cells constitutes the 
disordered topological conductor. It is connected to two 
perfect ballistic leads, one on each side. 

In order to determine the transfer matrix, the Hamil- 
tonian given by Eq. , which only connects sites within 
neighboring unit cells (cf. Fig. 14), can be written as 



^ ^ ^ ^n(H M;M _ iC M _i -|- H n?n c n 



H 



n,n+l^n+l ) i 



where 



(A.l) 



4 = (K,i) f (6 u ,i) f Ki) f ... (b u ,MY) 

is the vector of all electron creation operators in unit cell 



u, as indicated in Fig. 14 The matrices H UyW describe 



hopping from unit cell w to unit cell u. 

Starting with the Schrodinger equation one can find a 
recursive relation for the wavefunctions as 



^ u 



with the transfer matrix 



where 



A 

B 7 , 




A 



I 

B, 



H 



(A.2) 



(A.3) 



U,U+1 5 



: (H u ,ii-i) 1 {Epl — H u ,u) ■ 



M 



is a vector of 



the amplitudes on all sites in one unit cell and Ep our 
chosen Fermi energy. Using appropriate initial wavefunc- 
tions and Eq. ( A.2 ), we can calculate the wavefunction in 



every unit cell of our system. Our initial condition is an 
incoming right-moving wave in lead L such that there are 
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only transmitted waves in lead R and the initial and re- 
flected waves in lead L. Once we have the wavefunctions 
in both leads we can determine the transmission matrix. 
For computational convenience we calculate backwards 
through the system starting in lead R. 

The first step is therefore to determine the eigenmodes 
of lead R which are given by the eigenvalues A n and eigen- 
vectors \l/ n of the transfer matrix M w for w in lead R. 



For propagating modes 



ik n a 



and the eigenvalue is a 



root of unity, while otherwise the amplitude is increasing 
or decaying. Depending on the energy level we choose, 
we get p pairs of propagating modes (left and right mov- 
ing) with eigenvalues e ±%knCL (a is the length scale of our 
system and will be set to 1) and eigenvectors i &±k n and q 
pairs of exponentially increasing or decaying modes (re- 
sulting from imaginary momentum modes) with eigenval- 
ues e ±KnCL and eigenvectors i f?± Kn . The physical eigen- 



modes of lead R are the starting points for Eq. (A.2). 
Taking w to be a unit cell in lead R and u a unit cell in 
lead L, as illustrated in Fig. [l4j we can determine 2p + q 
physical wavefunctions in lead L using 



y u 



M 



u+1 



(A.4) 



and analogously for . These wavefunctions have to 
be a superposition of all possible eigenvectors of M in 
the leads and can be written as 



R ANk„ 



2q 

171=1 



, fc e"*" 



2p 
rn=l 

(A.5) 

and analogously for with N — w — u and k p + n = 
—k n: K q + n = — n n for brevity. OL q q > G C is the coeffi- 
cient of mode q m in lead L when mode q' n was propagated 
from lead R to lead L. 

Since we know that the eigenvalues of the transfer ma- 
trix are of the form e ±%k and e ±K , the momenta k m and 
can be calculated from the eigenvalues of M w . Eq. 



(A.5) forms a linear system of equations with a unique 
solution for the variables a a a > . In general all the co- 
efficients ot q q > will be nonzero, that includes the coef- 
ficients of the unphysical modes. Although we now es- 
sentially have corresponding incoming flux amplitudes in 
the left lead and outgoing flux amplitudes in the right 
lead, we cannot calculate the transmission coefficients 
from those because the picture is not physical. The wave- 
function in the left lead contains exponentially increas- 
ing contributions. Instead of choosing a single outgoing 
mode in the right lead (which by itself gives us unphys- 
ical contributions in the left lead) we take a superposi- 
tion of physical modes, i.e., propagating and exponen- 
tially decaying modes, in the right lead. The coefficients 
of this superposition have to be chosen such that (after 
propagating backwards through the strip) the unphysical 
modes in the left lead cancel each other. However this 
cancellation only works for a superposition of one propa- 
gating mode with all the exponentially decaying modes. 



lead L 


disorder 


lead R 
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FIG. 15: Considering a superposition of physical modes in 
lead R, one can cancel the unphysical modes in lead L by de- 
structive interference. The coefficients /3k m ,n have to be cho- 
sen such that, after backwards propagation, the exponentially 
increasing modes in lead L have zero weight. The coefficients 

Pkrr 



are given by Eq. (A.6). The coefficients 7q^,q^ of the 



resulting supe rpos ition (of only physical modes) in lead L are 
given by Eq. (|A.7). 



This procedure is illustrated in Fig. [15j We are not in- 
terested in the left-going modes in the right lead since 
we are looking at the case of an incoming wave from the 
left. To cancel the unphysical modes in the left lead we 
need to solve p systems of equations like (cf. Fig. 15) 



W 



(A.6) 

one for each k m corresponding to a mode propagating to 
the right. The /3fc m , n are the coefficients of the superpo- 
sition in lead R that we want to determine. The index 
k m denotes which propagating mode is part of the super- 
position. The index n shows which mode belongs to this 
coefficient, for the propagating mode and n = l...q 
for the q decaying modes. The ot q q i are calculated 
from Eq. (A.5) and are the coefficients of the unphysical 
modes that we want to cancel. C£ Kll k m is the coefficient 
of the (exponentially increasing) mode in the lead L 
when the (propagating) mode k m was propagated from 
lead R bac kwar ds through the system. The set of linear 
equations ( A.6[ ) is under determined and we can there- 
fore choose f3k m ,o = 1 without loss of generality. The 
coefficients Pk m ,o an d Pk m ,n gi ye us a superposition of 
physical modes in lead R which is caused by a superposi- 
tion of only physical modes in lead L. This is illustrated 
in Fig. [15] where the entire system is now made up out 
of only physical modes. In order to determine the trans- 
mission coefficients, we are interested in the coefficients 
7 of the propagating modes in lead L. Taking into ac- 
count the necessary superpositions, the coefficients for 
the right-moving modes are given by 



q 

7fcn ,fe m = 1 * ^k n ,k m + ^2 ft k ™,Ki ®k nj Ki 



(A.7) 
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lk n ,k m is the coefficient of mode k n in the left lead when 
the outgoing mode in the right lead is k m . The trans- 
mission coefficients can now be calculated by requiring 
probability current conservation. This leads to the fol- 
lowing system of equations: 



V k n — h Sf 



k=kn is the velocity of mode k n . The T nm 
are the transmission coefficients with \T^\ 2 the prob- 
ability that an electron coming in through mode k m is 
transmitted through the conductor into mode k n . 



rpLR 



\Tk R 



r = v, 



rpLR j 
... -Lpp / 



(A.8) 



with IV 



7k ni k m V\vkJ and 
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